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I.  INTRODUCTION 

Euler's  Tbeorem  Is  the  basis  for  a four  parameter  alternative  to 
Euler  matrix  representations  of  yaw,  roll,  and  pitch  motions  of  an 
unbending  missile.  The  quaternion  [1-3]  has  already  demonstrated 
superior  features  In  that: 

1)  It  needs  only  357.  of  the  multipliers  used  by  Euler  matrix 
computers. 

11)  Analog  and  digital  computations  are  more  accurate  than  Euler 
matrix  computations. 

Ill)  It  Is  preferred  to  Euler  matrix  representation  In  strap  down 
gyro  platforms  that  replace  expensive  mechanical  gimbals  with 
digital  computers. 

The  bordered  matrix  form  Is  emphasized  exclusively.  Transmitted 
quaternions  Interchange  lead  vectors  and  effect  speedier  hand  computa- 
tions. These  two  quaternion  Implementations  are  attributed  to  Ickes  [4] 

The  nth  root  of  a quaternion  is  an  instructive  precursor  to  small 
angle  recursions.  These  algebraic  recursions  are  regimented  by 
Chebyshev  polynomials. 

The  appendices  contain  four  glmbal  quaternions  for  strap  down 
computers  and  quaternion  representation  by  differential  equations. 


II.  SOME  PRELIMINARIES 


Slight  familiarity  with  Euler  matrices,  trace  of  a matrix,  and 
quaternions  will  be  presumed.  The  ordinary  matrix- vector  notation  and 
some  Gibbs  vector  (x,  y,  z)  algebra  will  be  employed.  The  transpose  of 
a matrix  or  vector  Is  denoted  by  a prime.  Unit  Gibbs  vectors  are 
Indispensable  and  are  denoted  by  a caret;  swift  Interpretations  of  the 
vectors  | and  | are  possible.  The  Euclidian  norm  ||||j  also  appears. 

Example  1 

If  I Is  a 3-vector  and  |||||  ■ s/PT>  then  | " l/IUH  . Other 
objects  such  as  the  lead  vector  and  transmitted  quaternions  will  be 
Introduced  as  needed. 


4 


ipmRiiiPi 


y,\MU^wv^y':  i 


III.  ROTATION  IN  THREE-SPACE 

Euler  matrices  produce  rotations  in  3- space  by  successive  rotations 
In  three  2-spaces,  the  xy-plane,^the  yz-plone,  and  the  zx-plane. 
Reference  vectors  are  x,  y,  and  z,  respectively,  the  first  symbols  In 
the  plane  designations.  The  plane  normals  are  z,  x,  and  y.  The  sign 
of  angles  is  determined  by  the  right  hand  screw  convention.  These 
rotations  produce  the  yaw,  roll,  and  pitch  angles,  respectively.  Euler's 
theorem  returns  a 3-space  rotation  to  its  basic  situation  in  2-8pace. 

Theorem  3.1  (Euler).  A 3-space  rotation  is  equivalent  to  a 2-space 
rotation  about  a fixed  vector.  The  2-space  rotation  has  the  fixed 
vector  for  its  normal. 

Sketch  of  Proof! 

a)  If  E ■ (®£j)  in  a 3 X 3 Euler  matrix  such  that  EE'  “I,  then 
there  exists  a fixed  3-vector  | such  that  E| 

b)  The  angle  of  rotation,  6,  is  given  by 


1 + 2 cos  e ■ ej^j^  + 622  + ©22  ^ tr  E (3.1) 

c)  Alternatively, 

cos  I - -I  yi  + tr  E , (3.2) 

sin  I - “ ,/tr(I  - E)  . (3.3) 

The  three  components  of  the  fixed  vector,  |,  and  the  angle,  9,  constitute 
a four  parameter  system.  Naturally,  there  exist  three  possible  classes 
of  representations: 

1)  Four  component  vectors, 

il)  Four  parameter  bordered  matrix, 

ill)  Four  parameter  pseudo- vectors. 

The  first  representation  will  be  considered.  Afterwards,  a return  to 
quaternions  in  bordered  matrix  form  is  in  order.  The  pseudo-vector  is 
sometimes  called  a hypercomplex  number.  It  is  Hamilton's  (1843)  original 
formulation  and  can  be  found  in  Bean  [3].  It  is  a matrix  decomposition 
of  the  quaternion  and  does  not  need  to  be  considered  in  this  computation 
context. 
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IV.  THE  CONSTRUCTIVE  EULER  THEOREM 

Theorem  3*1  can  be  reinforced  by  two  tasks: 

i)  Determine  two  vectors  In  the  2- space  whose  normal  is  the  fixed 
vector. 

11)  Construct  a four-component  lead  vector  that  deduces  a complete 
orthogonal  basis. 

As  Implied  In  the  statement,  a lead  vector  can  deduce  other  vectors  and 
with  them  form  an  orthogonal  (orthonormal)  bases. 

Lemma  4.1.  The  fixed  vector  | of  Theorem  3.1  satisfies  (E  - E')g  ■■  0, 
and  Its  components  are 


Moreover,  there  exists  a natural  vector  decomposition 
I « u - V , 

where 


Lemma  4.2.  The  vectors  (u  X v)  and  (u  x v)  X (u  - v)  are  in  the  plane 
whose  normal  Is  the  fixed  vector  |. 
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SCudlas  of  four  parametar  4-vactor  rapraaantatlon  can  commenca. 
Tha  norm  praaarvlng  mapping 


is  a formal  condensation  of  all  that  has  been  stated  so  far  into  the 
4-vactor  on  tha  right.  The  angle  Q is  implicit  in  the  Euler  matrix  E 
through  Theorem  3.1.  Norm  preservation  simply  means  that  the  4-sphere 
has  the  same  norm,  ||(||,  as  tha  3-8phare. 

/Hell  cos  A 

Lemma  4.3.  The  vector  I I is  a lead  vector  that  can  deduce 

\ 5 sin  0 / 

tha  other  three  members  of  an  orthogonal  basis.  The  orthogonal  basis, 
after  using  ( ■ u - v,  is 


II  u - v||  sin  6 
(u  - v)  cos  9 


.(u-  v)  X (uXv); 


If  c ^ u X v,  d ^ (u  X v)  X (u  - v) , then  the  orthonormal  basis  is 


\ I sin  9/  \ 5 cos  9/  \c/  \ d/ 

The  second  vector  is  obtained  by  replacing  9 in  the  first  vector  with 
9 *f  ivt/2).  No  more  angles  are  available,  and  the  top  components  of  the 
third  and  fourth  vector  are  zero.  Tha  scalar  triple  product  from  Gibbs' 
vector  algebra  ensures  orthogonality  of  the  vector  parts  of  the  third 
and  fourth  oasis  vectors. 

Lemnas  3.1,  4.1,  4.2,  and  4.3  comprise  the  Constructive  Euler  Theorem 
whose  corollary  is  an  orthogonal  basis  of  4«vectorB.  The  Constructive 
Bulor  Theorem  generates  the  need  for  4-vectorB,  but  computations  are 
slow.  Swift  computation  is  effected  by  quaternions  largely  because  cross 
products  are  present  in  quaternion  products. 
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V.  SIMPLE  QUATERNION  COMPUTATIONS 


The  quaternion  can  be  Introduced  through  the  nonstandard  Cayley- 
Kleln  fom 


(a  + jb  -c  + Jd 
c + jd  a - jb 

A nonunique  matrix  representation  of  complex  numbers  leads  to  the 
4x4  matrix 


(5.1) 


The  notation  Q(q)  means  the  "quaternion  array"  Q whose  first  column 
vector  Is  q.  The  standard  Cayley-Kleln  form  places  the  lead  vector  In  the 
the  first  row.  After  becoming  familiar  with  fundamentals,  a bordered 
4x4  matrix  will  be  used  exclusively;  this  usage  coincides  with  Ickes  [ 4]  . 


An  Intermediate  result  Is  needed. 

Lemma  5,1.  The  column  vectors  of  Q(q)  form  an  orthogonal  basis;  the 
column  vectors  of  Q(q)  form  an  orthonormal  basis.  Moreover,  q Is  a lead 
vector. 

Proof; 

a)  Q(q)  Q'(q)  - (q'q)  I (5.2) 

Is  obtained  by  direct  computation, 

b)  Q(q)  Q'(q)  - I.  (5.3) 


This  lemma  Immediately  replaces  many  geometrical  calculations  required 
by  the  orthogonal  basis  of  the  Constructive  Euler  Theorem  with  a format 
of  four  signed  and  unsigned  rearrangements  of  the  lead  vector  q that  are 
orthogonal. 


If  Indices  are  employed  for  the  lead  vector,  the  nonstandard  vector 
component  notation  (p^,  p, , P2,  p^)'  will  be  etqployed.  Oftentimes  Inter- 
est focusses  on  the  "vector  part"  (pj^,  P2*  P3) ' the  lead  vector  Is 
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shortened  to  (p^,  p) ' ; this  partitioned  vector  reduces  symbol  clutter. 
The  "real  part"  of  the  vector  Is  the  scalar  p^. 

Skew  quaternions  with  zero  real  part  are  Interesting  In  their  own 
right.  Direct  computation  gives 
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0 


(5.4) 


depending  on  whether  Gibbs'  vector  or  matrix  — vector  notations  are 
employed.  A computation  such  as  the  commutator  Identity 


requires  more  tools  for  effective  computation.  Two  digressions  are 
needed  to  consider  lead  vector  selection  rules  and  transmuted  quaternions. 

Finally,  the  bordered  matrix  form 

Pq  “P'  \ 

° (5.6) 

P PqI  + V 


exhibits  the  3x3  matrix,  (p^  I + P),  the  "kernel 


quaternion.  If  p^  ■ 0,  then  the  skew  matrix  F Is 


n 


Q|  1.  The  symmetric  part  of  Q| 
\p/  VP 


is  Q 


" Pr 


to  realize  that  ^ J automatically  constructs  the 
Interaction  of  real  part,  vector  part,  and  kernel 


" matrix  of  the 
the  kernel  of 

I.  It  Is  Important 
full  quaternion, 
should  be  apparent. 


VI.  LEAD  VECTOR  SELECTION  RULES 


Lead  vector  selection  rules  transform  matrix  multiplications  to 
simpler  matrix-vector  multiplications  and  thereby  reduce  calculation 
effort. 

Let  Dirac  brackets  (•)  denote  the  selection  of  the  first  column 
vector  from  a given  matrix;  l.e., 

^(p))  “ P . (6.1) 

The  following  rules  can  be  deduced  from  matrix  multiplication  rules: 


i) 

Distributivity  over  sums: 

{Q(p)  + Q(q))  - (q(p))  + (Q(q))  - p + 

CM 

• 

Sh 

ii) 

Dlstrlbutlvlty  over  products: 

<Q(P)  Q(q)>  = (Q(p)  <Q(q)>) 

(6.3) 

ill) 

Absorption  rule: 

<Q(p)  q)  “ Q(p)  q 

(6.4) 

iv) 

Idempotent  rule: 

{'(Q(p)))  ■ (p)  " p 

(6.5) 

Tlie  Idempotent  rule  Is  a slight  extension  of  the  absorption  rule  for  the 
case  where  Q iu  replaced  by  the  Identity  matrix.  In  the  dlstrlbutlvlty 
rule,  the  flow  of  Dirac  brackets  Is  from  outside  to  Inside  and  from  left 
to  right. 

The  Q operation  in 

Q(Q(p)  q)  - Q(P)  Q(q)  (6.6) 

Is  the  Inverse  of  the  absorption  rule. 

Many  of  these  rules  are  general  In  that  they  can  be  applied  to 
matrices  and  vectors  which  are  not  quaternions  and  lead  vectors. 
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VII.  TRANSMUTED  QUATERNIONS 


The  second  requirement  for  converting  matrix  multiplications  is  the 
Interchange  of  lead  vectors  on  the  left-hand  side  of  our  similarity 
equation 

Q(q)  P - Q<P)  q (7.1) 


to  the  preferred  lexicographical  order  p q on  the  right-hand  side. 

If  P'  ^ (••  b,  c,  d)  and  q*  ^ (a,  p,  7,  &) , 
then  form 


Q(p)  q 


/* 

b 

c 

Vd 


-b  -c  -d\ 

t\ 

/aa  - pb  - yc  - 6d\ 

a d -c  1 

pa  + ab  - 6c  + yd 

-d  a b 1 

ya  + 6b  + ac  - pd 

c -b  a/ 

w 

' 1 

^6a  - yb  + pc  + ad  / 

(7.2) 


Rewrite  the  rightmost  vector  so  that  Greek  letters  are  in  a matrix  and 
Latin  letters  are  in  the  column  vector  to  obtain 


Q(P)  q 


-y 

-6\ 

/‘\ 

e 

a 

-6 

y 

b 

y 

5 

a 

-P 

c 

V 

-r 

P 

aj 

'^yq)  p 


(7.3) 


If  one  suppresses  the  distinction  between  Greek  and  Latin  letters,  the 
array  is  obtained  from  Q by  transposing  only  the  kernel  matrix, 

Q^(q)  is  the  transmuted  quaternion  and  one  can  verify  that  q is  a lead 

vector. 

The  matrix-vector  result  can  be  summarized. 

Lemma  7.1. 

Q(q)  p - Q7(p)  q 

[Q7(p)]  [Q7(p)]'  - (P'P)  I (7.4) 

tr  Q^(p)  « tr  Q(p) 
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It  the  problem  is  interchanging  p and  q in  p'  Q(q)»  one  then  has  the  dual 


problem  that  leads 

to 

the  second 

transmuted 

quaternion.  Form 

/“ 

-7 

-8^ 

/oa  + 

pb  + 7C  + 6d\ 

a 

8 

-7 

f 

oib  - 

pa  + yd  - 6c 

p'  Q(q)  - (a 

b 

c d) 

-8 

a 

P 

■ 

0!C  • 

1 

pd  - ya  + 8b 

\6 

7 

-P 

al 

\ad  + 

pc  - yb  - 8a  y 

(7.5) 

Rewrite  the  rightmost  vector 

SO 

that  Greek 

letters  are  in 

the  row  vector 

and  Latin  letters  are 

In  the 

matrix  to 

obtain 

/* 

b 

c 

A 

b 

-a 

-d 

c' 

p'  Q(q)  - (a 

y 6) 

q'  Qi(p)  (7.6) 

c 

d 

-a 

* 

u 

-c 

b 

wA  / 

1£  one  suppresses  the  Greek  and  Latin  lettersi  the  array  is  obtained 

from  Q by  multiplying  the  second,  third,  and  fourth  column  vectors  by 
-1,  Moreover,  p is  the  lead  vector  of  Q,(p),  which  is  the  second  trans- 
muted quaternion. 

This  matrix-vector  result  can  be  summarized. 

Lemma  7.2. 

q'  Q(p)  • p*  Qj^(q) 

[q^Cp)]  [Qi(p)]'  - (p'p)  I . (7.7) 

tr  Qj^(p)  - -0,5  tr  Q(p) 

Ickes  [4]  employed  in  place  of  our  notation  and  called 
the  transmuted  quaternion.  This  notion  is  generalized  In  Appendix  A. 
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VIM.  FROM  QUATERNIONS  TO  EULER  MATRICES 

This  taction  reinforcet  pravlout  mathodt  to  achieve  firmer  founda- 
tions, The  first  task  Is  to  present  methods  from  behind  the  scenes  of 
Lemmas  3,1  and  4.1, 

Lemma  8.1.  Given  any  3-vector,  k.  It's  annlhllator K,  a 3-  x 3-skewmatrlx, 
can  be  constructed  as  follows 


Moreover,  If  g Is  another  3-vector,  then 
Kg  - -Gk  - g X k, 

K'  • K, 

G'  - -G, 
tr  - -2k'k, 

+ <k'k)  K - 0. 

Lemma  8.2.  If  h ■ Kg  and  K'  ■ -K,  then  H - KG  - GK  • gk*  - kg'. 

Proof  Is  by  direct  calculation, 

A frequent  problem  Is  that  K Is  known  but  g Is  unknown  and  subject 
to  K|  " 0,  The  solution  Is  | > k up  to  an  arbitrary  non-zero  scalar 
factor. 

Note  that  the  trace  operation,  trA,  as  the  sum  of  the  diagonal  ele- 
ments of  the  matrix  A. 

The  second  task^ls  to  derive  a more  specific  form  of  the  transposed 
lead  vector  (cos  d,  |'  sin  0). 

The  third  task  Is  to  derive  the  Euler  matrix  from  a quaternion. 

The  second  and  third  problems  are  solved  together. 


a unit  lead  vector,  then  the  similarity  transformation 


a unit  lead  vector. 


Pq  + p'p  " 1 

Induces  the  Euler  transformation 


[pp'  + (PqI  - r 4 E r 


Pq  - ±|cos 


p ■ the  fixed  vector  of  E. 


Proof; 


i)  Q 


(X)  ■ Q 


"p  po^  “ ^ p po^  - ^ r 


0 pp'  + (PqI  - P) 


2 U 


iv)  tr^pp'  + (PqI  - P)^)  - 3 Pq^  - p'p  - 1 + 2 cos  6 


implies  that 


[po^  - ("’•  I)  j [pq  ■ -I"”  (l)l 

v)  [^Q'(p)  QyCp)]'  [q'(p)  QyCP)]  ■ I implies  that  E'E  - I 

Vi)  (E|  - I and  EE'  - I)  imply  that  (E  - E')6.-  0 

vii)  (E  - E')|  - -4  PqP|  "0,  p^  ^ 0,  implies  i - p up  to  a 
multiplicative  scalar  factor; 


V-T"'  ■ ''■•Tif  VH'  •’ 


/Po\  . /“•(• 

\P/  \t>  "in 


vlll) 


Step  vlll  displays  an  amended,  principal  part  form  for  a lead  vector 
that  originally  contained  the  full  angle.  Note  that  the  ± sign  on 

In  Step  Iv  has  been  arbitrarily  made  positive  In  Step  vlll.  This  sign 
decision  must  be  made  external  to  the  quaternion  system. 

Some  straightforward  calculations  can  be  executed.  The  first  cal- 
culation Is  the  explicit  form  of  the  Euler  matrix,  E,  in  terms  of  the 
lead  vector  (p^,  p')  components;  namely 


y 

Lz  J 


2 2 
1 - 2p2^  - 2p3'‘ 

2P0P3  + 2piP2 

2P1P3  - 2PoP2 


2PiP2  - 2pjjP3 
1 - 2p^^  - 


2poP2  + Zpjpj 

2P2P3  - 2PqPj^ 

2 2 
1 - 2p^^  - 2p2  j 


(8.4) 


The  skew  matrix  Is 


m 

0 

-*V3 

‘PoPa  ’ 

0 P3 

“P2' 

-4PoP2 

0 

‘'oPl 

-4PoPi 

0 

0 

0. 

1 

■ 

-P3  0 

P2  -Pi 

Pi 

0 

(8.5) 


The  second  and  last  calculation  specifies  a product  of  quaternions 


<•)■<:)  <?)<:)  ■ 


(8.6) 


where  u,  v,  and  w are  mutually  orthogonal  vector  parts  to  greatly  reduce 
calculations.  Employing  x,  y,  z for  the  first,  second,  and  third  column 
vectors  of  the  Identity  matrix,  one  obtains  the  specific  form 


/ 

e ^ 

/ 

^ 

/ 

CM 

/ ^3  \ 

cos 

2 

- Q 

cos 

2 ' 

Q 

' cos 

2 

cos  Y 

p sin 

'V 

^ sin 

'a 

^ sin 

^9  .in 

(8.7) 
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The  strange  unit  vector  and  subscript  association  derives  from  our 
canonical  Euler  transformation  sequence;  i.e,,  rotation  first  in  the 
xy-plane  whose  normal  is  z,  rotation  in  the  y*-plane^ whose  normal  is 
X,  and  rotation  in  the  zx-plane  whose  normal  is  the  y.  Moreover,  the 
leftmost  quaternion  performs  the  first  operation  in  the  similarity  trans 
formation.  Direct  calculation  leads  to 


Pq 


- cos(f)-  (cos  (cos  (cos 

+ (sin  t)  («^"  t) 


(8.8) 


A 

P 


- (esc  f ) 


^08 

2 

. "2 
sin  -j- 

cos  Y + 

sin  y 

cos 

COS 

^^1 

2 

^2 
cos  -Y 

sin  y - 

, t^l 

sin  y 

sin 

is  in 

^^1 

2 

^2 

cos  -y 

^3 

cos  y - 

>^1 

cos  y 

sin 

^^2 

2 

^2 

T 

!l2 

2 


sin 


cos 


h 


sin  -«r 


2/ 


(8.9) 


It  should  be  emphasized  here  that  the  quaternion  angles  ^2»  s’f® 

the  negative  Euler  angles;  this  anomaly  arises  because  of  our  choice  of 
the  first  column  vector  as  the  lead  vector. 


IX.  SMALL  ANGLE  APPROXIMATIONS 

The  email  angle  approximations 

sin  ^ A |ji  I 

} for  I4 1 ;S  0.1  rad 
cos  ^ ill  1 1 


la  a mechanism  for  calculating  the  Euler  matrix,  E,  from  the  following 
Identity.  If 


K 


then 


i 0 

“•"3 

r2 


-Pi2 


^il 

0 


\ 


/ 


» 


(9.1) 


The  exponential  matrix  Is  also  the  culmination  of  the  differential  equa- 
tion formulation  of  rotations  (Appendix  C) . Although  differential  equa- 
tions can  be  solved  recursively,  one  would  rather  emphasize  direct  alge- 
braic methods. 


The  next  task  Is  to  compare  small  angle  second-order  effects  on 
Euler  matrices  and  quaternions.  Consider  the  quaternion  similarity 
transformation  again.  Quaternion  calculations  specifically  lead  to 
t|>e  Euler  matrix 

2 

pp'  + (PqI  - P)^  » I - (sin  e)P  + ^sln  |)  (pp'  + - I)  . (9.2) 

Some  simplified  calculations  are  In  order.  If 

- (1,  0,  0) 

and 

P'  ^ (0,  I,  1)  , 


then 


pp'  + p^  - I - -2  dlag  (p)  . 
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If  the  prevloui  Identity  and  the  email  angle  approximation  are  invoked, 
one  obtains 

pp»  + (PqI  - P)^  - I - et  - I diag  (p)  . (9.3) 

However,  if 

Cl)  ■ 

then 

2 

pp'  + (PqI  - P)^  - I - ep  - (I)  T , (9.4) 

T - diag(l,  -1,  -1), 

and  matrix  norm  comparisons  of  the  second-order  matrices  can  be  made; 
these  are 


(9.5) 


(9.6) 

One  can  conclude  that  the  quaternion  small  angle  approximation 
should  be  slightly  better  than  the  Euler  matrix  small  angle  approxima- 
tion if  the  goal  is  to  obtain  the  smallest  second-order  terms.  Such  a 
goal  is  implicit  in  the  exponential  matrix,  Equation  (9.1). 


II I diag  (p)  II  - (9^  - 0.7078^ 

III  T II  • ^ m 0.4330^ 


“fll-lnduc^i  r«ur.lon.,  it  l, 
of  2.  “ "l“«t»™lon.  Thl.  provide,  eo 

c^«lo°o  th*  l«r*. 

Hiree  mathematical  tools  must  be  used.  These  are  as  follows: 
s)  The  Cay  ley- Hamilton  Theorem* 

b)  The  minimal  polynomial  aquation  of  a quaternion. 

c)  The  modified  Chabyshev  polynomials  S (x)  and  C (x). 

n n' 


Cayley-Hamllton  - The  polynomial  equation  of 

[q^  - 2Pq  Q + (pq  + P'p)lj^  - 0 
and  the  minimal  polynomial  equation  for  a unit  quaternion  is 
" 2p„  Q - I 


is 


(10.1) 


(10.2) 


00.1.1  S_^(x)  with  thl.  .object  through  power,  of  . quaternlL  lo  the 
Identity 


- Vi«Po)  <(j)  - s„.a(V 

n - 2.  3.  4,  5,  . . . 

The  right-hand  side  of  the  equation  can  be  written 
^/o\  ^ ^O^n-l^^Po)  " S„.2(2Po)' 

Vn.l<2Po)  / 


(10.3) 


as 


(10.4) 
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ppi|>ilt»vy»»*H»t'<'.w%''"»<»wii'j>^'i> 


i 


or 


\p/  V S^.^(2Pq)/ 

The  last -equation  exhibits  the  remaining  Chebyshev  polynomial  C^(x). 

Lanczos  [5]  presents  a very  readable  Introduction  to  Chebyshev 
polynomials. 

Chebyshev  polynomial  Identities, 


C (2  cos  d)  " 2 cos  nd 
n'  ' 


and 


S , (2  cos  d) 
n-1 


(sin  ng) 
(sin  6) 


(10.6) 


(10.7) 


strongly  suggest  the  form  of  the  nth  root  of  a quaternion. 

Lemma  10.2.  If  n Is  an  Integer  and  n ^ 2,  then  the  nth  root  of 
/cos  0 \ / cos  (0/n)' 

Q(  )l8  Q( 

\p  sin  0 / \p  sln(0/n) 


Proof! 


0\ 

cos  — 1 

\ 

1 A 

1 

0.5  C (2 

1 

cos  0 

nV 

n/ 

1 

- Q 

ip  (sin  S 

.12  cos 

i) 

- Q 

ip  sin  0 

1 

\^\  n;  n- 

1\ 

V 

(10.8) 


In  the  similarity  transformation  context,  the  full  angle  0 must  be 
replaced  with  0/2, 

Computation  complexity  Increases  idien  cO8(0/n)  and  sln(0/n)  must 
be  computed  from  cos  0 and  sin  0.  Trigonometric  computations  first 
yield 
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cos 


■ 0.5  + cos  Q 


and 


cos^^^  “0.5 


1 + cos  ^ 


and,  in  general,  one  obtains  the  recursion 


cos 


«nri-l 


lu  * 0,  1,  2,  3,  ... 


(10.9) 


Finally,  trigonometric  computations  yield 


m “ 0,  1,  2,  3,  ...  . (10.10) 

Unfortunately,  Iterated  square  roots  In  both  computations  require 
multiple  precision  computations. 

The  nth  root  of  a quaternion  demonstrates  that  the  unit  quaternion 
property  must  be  maintained  in  a recursion.  Limitations  Inherent  In 
co8(€/n)  and  sln(d/n)  calculations  and  In  the  fact  that  n must  be  a 
power  of  two  can  be  avoided  by  the  Fade  approxlmant  to  be  presented 
In  the  next  section. 

With  Increasing  n,  1/n  tends  to  zero.  If  continued  further  to 
n ■ -1,  then  the  Inverse  problem  reoccurs  as  presented  In  Appendix  B. 
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XI.  SMALL  ANGLE  RECURSIONS 


Xhe  small  angle  approximations  are  restarted  with  the  norm  pre* 
serving  Fad£  approxlmant 


J(0/2n)  * 4n_+Ji  -^  < 0 1 rad  - -1 

* 4n  - je  • 2n  “ J ^ • 


which  implies  that 


2 2 
. 16n  ■ e 

"22 

16n^  + 0 


0 . 8ne 

’TTT—i 

Ion  + 0 


(11.1) 


(11.2) 


(11.3) 


If  0 ■ 1 rad  ■ 53  deg  and  n 
In  this  case, 


5,  then  |d/2n|  is  at  the  equality  boundary. 


'(«)* 

a* 


40 

401  * 


and  the  computer  binary  words  must  contain  at  least  9 bits, 

/ cos (0/2)  \ 

The  nth  root  of  Qi  ] is  approximated  by 

\p  sln(0/2)/ 


2 2 
16n^  - e 


2 ^ g2  8n0  p 


The  finite  matrix-vector  iteration  format  is 


Entering  these  values  Intu  the  Euler  matrix  and  using  the  Identity 

rl  +al«‘  0.25  C^(2x)  + (1  - x^  ,(2x)  - 1 (11.8) 

n n n n“  i 

yields  a slight  simplification  In 

pp'  + (PqI  - P)^  - I - 27^0^  P + a^(pp'  + P^  - I)  . (11.9) 

A slight  hazard  occurs  when  ^^(2?^)  = 0;  this  results  in  much  compu- 
tation to  produce  a simple  result, 

These  end  results  suggest  that  dedicated  digital  computer  Itera- 
tions may  be  emulated  by  the  Iterative  processes  of  C^(2Pq)  and 
_ 

1 - Po  » thereby,  matrix-vector  Iterations  and  normalizations 

are  minimized.  The  sole  matrix-vector  calculation  appears  In  the 
quaternion  to  Euler  matrix  conversion  near  the  end. 

The  minimum  multiplication  recursion  for  the  modified  Chebyshev 
polynomials  [5]  is 

So(x)  - 1 
S^(x)  - X 

S^(x)  - X S^_^(x)  - » (11.10) 

followed  by  additions  In 
Cq(x)  - 2 

Cj^(x)  - X 

^n^*^  “ ®n^*^  ■ ®n-2<’'>  • (11.11) 
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END 


Figure  1.  Single  quaternion  small  angle 
recursion  and  conversion  to  an  Euler 
matrix. 


XII.  ITERATION  OF  PRODUCT  QUATERNIONS 

Suppose  the  quaternions  producing  yaw,  roll,  and  pitch  motions 
In  the  similarity  transformation  are  symbolized  by  A,  B,  and  C, 
respectively,  and  that  the  combined  transformation  Is  the  product  ABC. 


In  Section  XI,  (ABC)' 
yielded 


was  the  small  angle  approximation  and  Iteration 


(ABC) 


The  exponential  matrix 


/O  0'\ 

» 1 * 

\0  K / 


of  Section  IX  In  the  context 


of  a triple  matrix  product  uses  an  extremely  large  Integer  n such  that 

K in  [I  + (K/n)]"^  retains  only  first-order  terms  and  the  iteration 
converges  to  the  Euler  matrix.  If  n Is  not  large,  then  a fundamental 
Inequation 


(ABC) 


1/n  ^ ^1/n  gl/n  ^1/n 


(12.1) 


appears  behind  the  scenes  In  the  quaternion  aspect. 

The  third  possibility  Is  motivated  by  the  desirable  equality 


" ABC 


(12.2) 


This  iteration  requires  three  parallel  Iterations  and  a final  matrix 
product.  This  results  In  a maximal  computer,  but  would  give  the  safest 
computation  for  the  finite  arithmetic  available  on  a dedicated  digital 
computer. 

A flow  chart  summary  of  the  primary  results  obtained  Is  presented 
in  Figure  2.  The  starting  equation  Is  the  quaternion  similarity 
transformation.  Ascending  or  descending  vertical  arrows  represent 
major  transitions  whose  meanings  are  apparent,  while  horizontal  arrows 
represent  rather  minor  transitions.  ENDS  is  a theoretical  ending, 
whereas  END4  and  END5  represent  computer  configurations.  A solidus  (/) 
through  a directed  line  denotes  an  unrecommended  transition.  END  numbers 
are  numerically  ordered  according  to  the  occurrence  of  results  In  this 
report.  Unrecommended  exits  are  assigned  the  last  numbers,  l.e.,  6, 

7,  and  8. 


END4  and  EHD5  ara  posslbla  nulcipla  axlti  for  this  study.  END4 
laads  to  axiatlng  dadieated  quatarnlon  e«nputars.  If  specifications 
should  tighten  In  the  future,  then  ENDS  Is  available  for  alternative 
consideration. 
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Xill.  CONCLUSIONS 

Datalled  quaternion  applications  to  three-dlmenelonal  rotations 
of  a rigid  body  have  been  presented  In  a bordered  matrix  format.  There 
Is  an  initial  abstruseness  about  quaternions  that  eventually  becomes 
more  concrete.  Recursions  can  be  generated  algebraically  after  Intro- 
ducing the  nth  root  of  a quaternion.  A flow  chart  for  the  quaternion 
similarity  transformation  study  summarizes  procedures  and  results. 

Two  appendices  contained  in  this  report  are  concerned  with  strap 
down  digital  computers  modelled  by  four  glmbel  quaternions  and  quater- 
nion representation  by  a matrix  differential  equation.  Numerical 
Integration  of  differential  equations  Is  understood  by  experts;  there- 
fore, because  of  this  and  emphasis  on  simpler  algebraic  recursions, 
no  differential  equation  (DE)  recursions  were  presented.  Quaternion 
fundamentals  range  over  several  mathematical  subjects  even  though 
angular  acceleration  of  a rigid  body  has  not  been  Included. 

The  reference  list  samples  the  extensive  quaternion  literature. 
Branets  [6]  has  been  oientloned  because  of  the  large  list  of  European 
papers,  books,  and  static  and  dynamic  applications. 

This  report  is  oriented  coward  analog  and  digital  computation  and 
fosters  fundamentals  which  eventually  produce  the  following: 

1)  Improve  computer  hardware  acceptance  testing. 

11)  Provide  better  formulation  of  computer  hardware  requirements. 

Ill)  Provide  minimal  software  emulation  of  hardware  realizations. 
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Appmdix  A.  GENERAL  QUATERNION  TRANSPOSES 

The  two  transmuted  quaternions  can  be  Imbedded  In  a set  of  eight 
general  quaternion  transposes.  Also,  there  exists  a natural  subscript 
notation  with  simple  binary  number  operational  rules.  Motivation 
arises  from  the  fact  that  transposes  of  quaternions  are  equivalent  to 
transformations  with  a diagonal  matrix  T. 

A-l.  If  T ■ dlag(l,  -1,  -1,  -1)  and  the  prime  denotes  matrix 
transpose,  then 

Q'(q)  - Q<Tq) 

and 


In  general,  T can  also  be  a left  multiplier  or  a right  multiplier  of 
Q as  In 


QjCq)  - T°Q(T°q)  t1  - Q(q)  T 

The  usual  binary  representation  of  1,  [1]2  " 001,  Immediately  locates  T 

In  a position  accordingly  because  1 Implies  the  presence  of  T and  0 
Implies  the  absence  of  T In  a particular  bit  position. 

It  Is  natural  to  consider  the  subscript  as  an  operation;  namely, 
the  generalized  quaternion  transpose 

Q^(p)  - [QCp)]^ 

The  next  task  Is  to  determine  k In 

If  1,  J are  given  from  the  set  of  Indices  0,  1,  2,  3, 4,  5,  6,  7. 


WPIPP 


' Tn''T’’r  I “ 


Lemma  A- 2.  If 

[l]2  - a^a^aQ  . 

[JJ2  " ^£^1^0  ’ 

[k]2  - <2'^iKq  , 

and  V Is  the  term-by-term  sum-modulo-2  of  a pair  of  binary  words,  then 

1)  © 32)(«i  ® ^iX^o  © V " ^ 

11)  k ■ 2(2K2  + K^)  + Kq 

gives  the  decimal  number  k in  [Qj^(p)  ■ Qj^(p). 

Table  A-1  displays  k for  given  1 and  J,  Here,  it  is  desirable  to 
shorten  the  precise  [ui,  V [J]2j  ^ typical  dogleg  path  for 

determination  of  k ■ 4VJ  is  shown  by  a pair  of  heavy  lines.  Some  simple 
arithmetic  representations  of  IVJ  are  also  noted  in  Table  A-1.  It  is 
natural  to  consider  the  fundamental  subset  of  generalized  transposed 
quaternion  as  Q^Cp),  Q^(p),  and  Qj(p). 

Lemma  A- 3.  The  complete  set  of  transposed  quaternions  can  be  determined 
from  the  fundamental  subset  Q^,  Q^,  and  Q^.  The  complete  set  can  be 
conducted  as  follows: 

i)  Q^»  Qj^  > 

li)  , [QjJ  yield  Q^,  Q^,  and  Q^,  respectively; 


ill)  Q2  " ^0  " ^ Implicit  in  the  construction  of  1) 

and  11). 

The  quaternion  meanings  of  Qy(p)  and  Qj^(p)  have  already  been 

presented  in  a lead  vector  interchange  context.  The  meaning  of  Q^(p) 
is  that  it  is  the  row  transformation  counterpart  of  Q^(p). 


i-  |-j  liYii'iiil't 
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The  subscript-?  operation  Is  employed  most  frequently.  Its  transpose- 
like nature  Is  now  demonstrated. 

Lemma  A- 4. 

[q(p)  Q(q)]7  - QyCq)  QyCp)  - Q^Cq)  QsCp)  • 

Proof; 

[Q(p)  Q(q)]7  - T[Q(p)  Q(q)]'T 

- [T  Q(Tq)]  [Q(Tp)T],  (T^  - I) 

■ Qg(q)  Q3(p) 

- (T  Q(Tq)T)  (T  Q(Tp)T) 

" QyCq)  QyCp) 

Only  subscripts  2,  3,  6,  and  7 have  middle  ONEs  In  their  binary  represen- 
tations and  these  subscript  operations  "transpose"  matrices  as  In  Lemma  A-4. 
Further  Investigation  can  proceed. 

Lemma  A-5. 

tQ(p)  Q(q)]2  " 

lQ(p)  Q(q)]3  - 

[Q(p)  Q(q)l6  - Q^Cq)  Q2(p) 

As  before,  these  factorizations  are  not  unique. 
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Appendix  B.  FOUR  QiMBAL  QUATERNIONS 


Quaternion  inverse  equation  variants  lead  to  two  error  quaternions 
applicable  to  the  remote  control  of  two  rotating  rigid  bodies.  A further 
consequence  is  the  four  glmbal  quaternion,  a mathematical  model  for  a 
strap  down  digital  computer  that  replaces  expensive  mechanical  gimbals. 


The  motivating  quaternion  equation  is 


Q'l 


( 


cos  6 


p sin  e 


(B-1) 


Replace  the  first  lead  vector  with  another  lead  vector  and  form  the 
quaternion  difference 


0/  cos  \|r  \ / cos  6 \ 

^ \Jq  sin  i ) \ sin  9/ 


where  (cos  9,  p sin  9)'  Is  the  reference  lead  vector. 


Lemma  B-1. 


/I  \ / cos  ijr  \ / cos  9 \ 

W/  ^ \q  8in  Ik/  V sin  9/ 


i 


1 - cos  ijf  cos  9 - (q  • p)  sin  \k  sin  9 
q sin  i|r  cos  9 - p cos  >|r  sin  9 + (p  X q)  sin  t sin  9i 


"(s 


cos  9 - cos 


sin  9 - q sin  \|r 


) 4:°' ' ) 

j \psin9/ 


(B-2) 


Understanding  can  be  furthered  by  making  a side  calculation  for  the  special 
case  q " p.  Differences  of  angle  will  appear  in  such  a way  that  all  lead 
vector  components  aro  zeroed  when  the  two  rigid  bodies  are  aligned  and 
have  Identical  anglej.  This  result  reinforces  Ickes'  [4]  practical  goal 
of  applying  quaternions  to  aligning  two  rotating  systems. 


Replacing  one  lead  vector  with  an  orthogonal  lead  vector  In  the 
left-hand  side  of  Equation  (B-1)  leads  to  a complementary  Identity, 
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cos  d 


)cr:)-c) 


sin  9 


Transposing  both  sides  yields  the  following  result. 


(B-3) 


Lemma  B-2. 


sin  6 


H 


cos  9 


p cos  9/  \p  sin  9j 


(B-4) 


and  the  quaternion  product  commutes.  The  second  error  quaternion  comple- 
ments the  quaternion  previously  derived. 


Lemma  B-3t 


1)  - i‘"  * ) i 

\p/  \q  cos  t/  Vp  sin 


cos  9 


- sin  >lr  cos  0 + (q  • p)  cos  i|f  sin  9 
-q  cos  t cos  9 + p(l  - sin  sin  0)  + (p  X q)  cos  ^ sin  0 


(sin  0 - sin  t 
p cos  0 - q cos 


cos  0 


p sin  \lr 


(B-5) 


where  (cos  0,  p sin  0) ' Is  the  reference  lead  vector.  This  again  leads 
to  an  all  zero  lead  vector  when  the  two  rigid  bodies  are  aligned. 

An  additional  application  of  Lemma  B->2  Is  the  nulling  of  the  earth's 
spin  angle,  \|r  cot,  along  the  spin  axis,  r,  for  a long  range  missile's 


geonavlgatlon  system.  If  Q| 


O' 


In  the  original  similarity  transforma- 


tion, Is  replaced  by  a counterpart  of  Equation  (B-3),  one  obtains  the 
aslmllarlty  transformation 


0 ) 
cos  J ] 

sin^ 


/ sin  \ / cos  lie  \ / f \ P\ 

cos  J % sin  J \ sin  |/  " \p/ 


<B-6) 
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L«ad  vaetor  salaction  and  lead  vector  Interchange  on  the  aelmllarlty 
tranaformation  yield  the  dailred  reault. 


Leiana  B-4.  The  aelmllarlty  tranaformation  yields  the  four-glnibal  quater- 
nion model 


(B-7) 


% ^ 
which  nulls  the  spin  angle,  i|r,  along  the  spin  axis,  r.  Three  of  the  four 

glBd>als  are  resident  In  the  half-angle  quaternions,  whereas  the  fourth 


ginibal  la  represented  by 


The  extraordinary  event  here  Is 


that  the  angle  t appears  on  the  left-hand  side  but  disappears  on  the 
right-hand  side  — a hidden  variable  situation.  From  this  point  on, 
p Is  the  same  as  given  In  the  text. 


Another  scenario  for  four-glmbal  quaternions  is  the  spin  stabilized 
missile.  Missile  spin  produces  desirable  aerodynamic  stability  but 
undesirably  spins  radar  homing  data*  Four-glmbal  quaternions  "de-spln" 

the  data. 


The  simple  flowchart  for  error  and  four  glnd>al  quaternions  Is  shown 
In  Figure  B-1. 
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(END1) 


(END  2)  (ENDS) 


Figure  B-1.  Flow  chart  for  investigation  of  two  error 
quaternions  and  the  four  gimbal  quaternions. 
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Appandix  C.  QUATERNIONS  AND  DIFFERENTIAL  EQUATIONS 

The  matrix  DB  aspect  of  four  gimbal  quaternions  is  presented.  A 
simple  computation  is  considered  first. 


Lenana  C~I.  If 
i(t)  - X(t) 


0 


and 


X(0) 


idiere  r is  a constant  3-vector,  u « 11^11,  X(t)  is  a 4 X 4 matrix  function 
of  time,  0 j t < <»,  then 


1)  X(t) 


ii)  <X(t)> 


/ cos  J \ / cos  U)t  \ 

^p  sin  “/  sin  wt  / 

■i-M 


cos  U)t 


r sin  u)t, 


Proof: 


1)  1^  r^^  + r'r  ■ r*r 


2)  Q‘ 


0”“''  • 


3)  X + w^X  - 0 


cos  ut  + ^ j » 

('  COB  J \ / cos  «t  \ 

f ainf/'^V.ln.ty  ’ 


4)  X(t)  - X(0)  I I cos  ut  -I*  Ql 

e 

/ cos 

5)  X(t)  - Q 
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6) 

7) 


<X<t)>  - Qy 


End  of  proof. 


It  is  apparent  that  the  matrix  and  vector  portion  of  the  four-glmbal 
equation  with  >|r  (tH:  have  been  represented  by  matrix  differential  equa- 
tions. It  suggests  that  the  two  leftmost  matrices  be  computed  next 
as  the  slightly  more  difficult  case. 


Lemma  C-2.  If 

Y(t)  - Y(t)  Q| 

and 

Y(0) 


C) 


■'( 


cos  J 
sin  fy 


where  r Is  a constant  3-vector,  u)  * Ill'll*  0 a*  t < »,  then, 

cos  u>t 


1)  Y(t)  - Q' 


6 

cos  J 


- 4 9 

p sin  j 


Vr  sin  wtj 


t / cos  I 

11)  -w  / YCt)  dT  + Q'(^  0 ) Q( 

o yp  sin  ^ 

9 


cos  T \ (—  sin  wt\ 

" ^ . 0 / / 
^p  sin  jy  yr  cos  ut / 


Executing  step  11)  of  Lemma  C-2  and  step  11)  of  Lemma  C-1  as 
parallel  computations  followed  by  matrix-vector  multiplication  yleldE 


sin  ut' 


r COP  ut. 


the  end  result  of  a matrix  dlfferemtlal  equation  representation  of  a 
four-glmbal  platform. 

It  Is  now  possible  to  branch  Into  numerical  Integration  according 
to  Barker  [7]  or  Into  analog  simulation  according  to  Mitchell  [8]. 
Critical  comments  are  In  order  for  Barker's  [7]  numerical  methodology. 
Barker  compares  the  self-starting  Adamo-Bashforth-2  method,  which  belongs 
to  numerical  transform  methods  peculiar  to  electrical  engineering,  with 
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